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ABSTRACT 


This thesis investigates the accuracy of a recently proposed passive bearings-only 
Target Motion Analysis (TMA) procedure. The primary method of analysis 1s to com- 
pare computer generated positions of a Target that 1s moving with a constant course and 
speed, with the procedurally derived estimated positions. 

A computer model was developed which simulated several possible interactions be- 
tween the Target and Own Ship. Estimated parameters of the Target track were com- 
puted using the procedure under analysis. These values were compared to the “true” 
values generated by the simulation. 

An analvsis of the TMA procedure as originally proposed showed that it failed to 
accurately estimate the target track parameters. However with some modifications, the 
accuracy improved significantly and it is felt that the procedure can accurately estimate 


target range (but not necessarily course and speed) for some target geometries. 
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I INTRODUCTION TO THE PASSIVE RANGING PROBLEM 


A. BRIEF HISTORY 

The importance of doing Target Motion Analysis (TMA)! by a passive bearings- 
only tracking technique has had the attention of many tacticians over the last 45 vears, 
mainly because the submarine’s objective is to remain covert (undetected) during the 
tracking maneuver. This guarantees the advantage of surprise if an attack is to be per- 
formed. 

In the early part of WWII, TMA by passive bearings-only methods began to be 
used, and the following drawn from Naval Tactical Decision Aids bv D. H. Wagner, 1989, 


provides a brief history of its development: 


Munn W Wil LY. F. C/ Linch developed the Linch Plot. which uses a pivotal 
relationship among bearings, bearing rate, and Target relative motion. Based on 
that. he developed a graphical method that was used until the late 1960's. It was 
abandoned because it was useful onlv for the short ranges contemplated by Linch 
at that time. 


During the 1950's, various human plot methods and nomographs were used for 
ieee ne of the oldest of these was the Strip Plot, later called Geographical Plot. 
This method has often been used to provide upper and lower bounds on Target 
range. based on bounds on speed. An additional much used nomographic device 
has been the Bearing Rate Slide Rule. An early TMA computer was the Position 
Keeper. which was a carrvover from WWII. and used as an aid to approach and 
attack surface ships. 


In 1953, F. N. Spiess published a graphical procedure requiring four bearings and 
a maneuver by Own Ship. This method is known as Spiess Plot. It remains in use 
in surface ship TMA. 


The 1954 command thesis of LT. J. F. Fagan derived a four-bearing TMA sol- 
ution, Which was a transcendental system of three equations with three unknowns. 
To reduce to computability by slide rule. he assumed that Own Ship motion during 
the first three bearings was approximately zero which was probably satisfactory for 
diesel operations. 


In 1957, LT. J. J. Ekelund devised one of the most famous TMA methods, now 
called Ekelund Ranging. 


Between 1967 and 1969. D. C. Bossard with a group of USN officers, made an 
analysis of the fundamentals of bearings-only ranging and the effect of Own Ship 


] Determination of Target range. course, and speed 


maneuvers on ranging accuracy. Through this work and trials at sea, they observed 
that ranging errors could be eliminated by judicious choice of the time for which the 
range was estimated. Also, by selecting some particular maneuvers, this best time 
could be controlled to be in the past, present, or future; then they applied this 
method to various forms of bearings-onlv ranging, especially to Ekelund and Spiess 
methods. The general name for these procedures is Time Correction. 


In August 1968, a series of exercises was conducted aboard the USS Pargo to 
compare the different methods in use at that time. The Passive Ranging Manual 
evolved in three volumes from the report of the analysis of these exercises and was 
published later as an NWP. This 1s a basic reference today for TMA. 


In the late 1960's the CHURN method appeared, which uses a large number of 
bearing observations and regression methods to do TMA. This method was devel- 
oped by General Dynamics’ Electric Boat and Librascope. 


In 1969, H. W. Headle, J. Di Russo, and i. Viessere developed the (ianua 
Adaptive 1MA Evaluator (MATE). This method begins with an input estimate of 
Target course, speed and range, and from there estimated bearings are computed 
and compared against observed bearings. Differences between these bearing values 
are plotted against time. If these differences are distributed roughly symmetrically 
as white noise about zero. while Own Ship changes course or speed, then the trial 
solution 1s a good one; otherwise the MATE operator adjusts the input Target pa- 
rameters. 


In the early 1970's, the application of Kalman filtering appeared for the first time 
in TMA. It was a done by J. S. Davis of NUSC, based’ on the Pin es thesis aimee 
J. Murphy at Northeastern University}. Vater versions were devcloscatoy tn. 


In the mid 1970's, a new TMA method known as FLIT was devised by ENS. L. 
Anderson which has seen considerable effective operational use. Its methodology 
remains classified. 


In the mid and later 1970's, ranging and tracking at long range on a sphere were 
developed by D. C. Bossard, J. B. Behrla, and LK. Graves ci Daniele ene m 
Associates. This was motivated by over-the-horizon targeting needs. 


In the late 1970's and during the 1980s, further developments were made by using 
Kalman filtering and also stochastic differential equations. The Maneuvering Target 
Statistical Tracker (MTST) is a method devised from this basis by W. H. Barker of 
Daniel H.Wagner, Associates. 


From 1985 to 1987 the Generical Statistical Tracker (GST) was developed at 
Daniel H!1. Wagner, Associates. This is a PC version of MTS1. 
B. PURPOSE OF THE THESIS 
In laré 1988, LCDR. P. kK. Peppe of the Uss ia lella (Ss 70l) prencscam. 


bearing-only Target Motion Analvsis (I M.A) procedure (here called the Bearings Ex- 


trapolation Procedure) for use in submarines [Ref. 1]. It is the purpose of this thesis to 


test the accuracy of this procedure against a nonmaneuvering Target. 


C. PROBLEM DESCRIPTION 

A submarine (Own Ship) is trying to determine the range, course, and speed of a 
ship (Target) that is moving with a fixed course and speed. The data collected by Own 
Ship consists of passive bearings. Since the bearing data is obtained by a single tracker, 
a unique tracking solution may only be obtained by a maneuvering Own Ship [Ref. 2]. 
A TMA method which uses a two-leg maneuver can fulfill this requirement. 

The two-leg TMA procedure used 1n this study is depicted in Figure 1. Own Ship 
begins by receiving bearings to the Target each 20 seconds for 4 minutes. In Figure 1, 
i, 1s the time at which bearing / is received. After time 2, the bearings stop and Own 
Ship begins a course change. At time 4, Own Ship is on the second leg and bearings 
resume. Bearings continue on the second leg, the last one being received at time 4). 
Estimates of Target course, speed and range will be derived from these 24 observed 
bearings. 

The central idea of the Bearings Extrapolation Procedure 1s that given three bearings 
observed at three times on a single leg, all future and past bearings can theoreticallv be 
determined. (This idea was known bv Spiess [Ref. 3] and mav predate him.) So, refer- 
ring to Figure 1, the estimated Target position at time é can be obtained by crossing the 
true bearing (beginning at point 4) with the extrapolated bearing from leg two (begin- 
ning at point 4’). 

The same as above can be done for all time points that have bearings between 4 and 
f,. Then regressing linearly on these estimated Target positions, the Target track pa- 
famevers Cal) be determined. 

Although this procedure should give exact answers, in practice large errors can be 
generated. The reason for these inaccuracies is that bearings to the Target cannot be 


measured exactly, so the extrapolated bearings will also be in error. 


target track 


fix ot (41 


F. 


/ bearing 
/ at (47 
® id 
< | il 
SY ee A” tat" 
own ship : vue 
N 


second leg aes 7 133" 
a \ 


extrapolated 


extrapolated 
legs 


“ VN 
own ship aoe 2 


first leg t6' ~ 
— t3' pe 





Figure 1. Graphical Representation of the Bearings Extrapolation Procedure. 


Exactly how bearings that contain errors are extrapolated, with as much accuracy 
as possible, forward and backward in time 1s part of what this thesis examines. The most 
straightforward method (and the first to be analyzed) is simply to use linear regression 
on the (B, , t) data, where B is the bearing at time t. Alternate methods include re- 
gressing the (8, , t) data by using quadratic or arctangent models. Each of these proce- 


dures were tried and accuracies assessed. 


The Spiess Graphical Method, Ekelund Ranging, and the CHURN Method are 


closely related to the procedure proposed by Peppe, and for that reason they are further 


explained in Appendix A. 


D. 


tJ 


a 


tJ 


ASSUMPTIONS FOR THE SIMULATION 


The computer simulation used in this study uses the following assumptions: 


The Target moves at a constant course and speed throughout the tracking maneu- 
er. 


Differences in depth for the Target and Own Ship are disregarded. 
Errors in position for Own Ship are disregarded. 


Errors in bearings to the Target are independent and identically distributed normal 
random variables with mean zero and standard deviation not larger than 1.5°. In 
addition, this standard deviation is constant for each simulation experiment. 


Bearings to the Target are obtained each 20 seconds. 


Initial course, speed and range for the Target and initial course and speed for Own 
Ship are as described in Chapter II. 


TERMINOLOGY 


Lead-leg. An Own Ship course such that Own Ship speed across the line of sound 
is in the same direction as the Target speed across the line of sound at the start of 
une leg, 


. Lag-leg. An Own Ship course such that Own Ship speed across the line of sound 


is in the opposite direction to the Target speed across the line of sound at the start 
Oretine leg. 


Leg length. The distance along either leg where bearings are received (does not 
include Own Ship turn). 


Raw bearing. An unsmoothed bearing to the Target. 


Faired bearing. A bearing obtained after smoothing of the raw bearings. 





Extrapolated positions. Positions of Own Ship obtained by extending the first leg 
forward in time or the second leg backward in time. 


Extrapolated bearings. Bearings obtained by extrapolating actual bearings forward 
in ume for the first leg or backward in time for the second leg. 


Lead-leg estimation. The estimated target track obtained bv linear regression of 
estimated Target positions. These estimated Target positions are obtained from the 
intersection of faired bearings from the lead-leg and extrapolated bearings from the 
lag-leg. 





Lag-leg estimation. The estimated target track obtained by linear regression of 
estimated Target positions. These estimated Target positions are obtained from the 
intersection of faired bearings from the lag-leg and extrapolated bearings from the 
lead-leg. 


10. Lead angle. Angle between the line of sound and the course of Own Ship when it 
is leading the Target. 


11. Lag angle. The angle between the line of sound and the course of Own Ship when 
itis lagging the Target. 


12. Angle on the bow. The angle between the line of sound and the course of the 
iV aRees. 
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Il. TESTING THE BEARING EXTRAPOLATION PROCEDURE 


A. INITIAL CONDITIONS 
In order to perform a thorough analysis of the Bearing Extrapolation Procedure, 
5,184 separate simulation experiments were conducted. Each experiment used one of the 


possible combinations of the following parameters: 


1. The initial Target positions were all due north from Own Ship (i. e., 000°) at one 
of three different initial distances: large (= 60,000 yds), medium (= 30,000 yds), 
and short (= 10,000 yds). 


feeoix courses for the Target were used: 030°, 060°, 090°, 120°, 150°, and 173°. 
emmlnree speeds for the Target: 25, 15, and 5 knots. 

4, Four courses for Own Ship in the lead-leg: 060°, 070°, 080°, and 090°. 

5. Four courses for Own Ship in the lag-leg: 270°, 280°, 290°, and 300°. 

6. Three speeds for Own Ship: 5, 10, and 15 knots. 

7. Two maneuvers by Own Ship: lead-lag and lag-lead. 

8. Four bearing error standard deviations: 0°, 0.5°, 1°, and 1.5”. 


It is important to mention that because of problem symmetry and the fact that both 
lead-lag and lag-lead maneuvers are simulated, there was no need for simulating courses 
for the Target from 180° to 360°. 

The randomness in the simulation was introduced only through errors in the re- 


ceived bearings. 


B. DESCRIPTION OF THE SIMULATION 

The simulation was written in Fortran 77 [Ref. 4] and executed on an IBM 30353 
mainframe. Appendix B has a complete listing of all variables of the simulation pro- 
gram. Appendix C is a listing of the simulation program code. 

For each combination of initial conditions, the simulation was repeated 100 times 
with different seeds for pseudorandom-generation of bearing error. For each of the re- 
petitions, the following steps were conducted: 

1. Target positions Were simulated every 20 seconds for the entire maneuver. 
2. For the first leg the following data was generated: 


CaO wes position every 20 seconds. 


b. Bearing to the Target everv 20 seconds (normal errors are added). 


c. Own Ship track extended bevond the turn. 
3. For the first leg the following data was then computed: 
a. Coefficients for linear and quadratic fits for bearings versus time. 


b. Faired bearings and extrapolated bearings using the results from linear or 
quadratic fit. 


4. Own Ship course was determined for the second leg; based on the required second 
leg. lead or lag angle. and the line of sound specified by the extrapolated bearing 
at the nuddile of the imaneuycmirn): 


5. Simulations and computations for the second leg were completed in the same 
fashion. 


6. Estimated positions of the Target in two coordinates (X and Y) were computed by 
deternuning the point of intersection of faired bearings with their corresponding 
extrapolated bearings. 


7. Linear and orthogonal regression 1s conducted on the estimated Target positions 
(in X versus Y coordinates), to determine the Target course. 


8. Range to the Target at the end of the maneuver was obtained by computing the 
distance between Own Ship’s actual position at the end of the second leg and the 
estimated Target position. 


9. Finally, Target speed was obtained by computing the distance between the esti- 
mated Target position at time 4s, and the estimated Target position at the end of the 
second leg (4). on the fitted Target track, and then dividing this distances@: 
ls) i a [y. 

The simulation of the bearing errors was done using the Linear Congruential 
Method for unirorm pseudorandom numbers, combined with the Box and Muller 
Method to produce normal vaniates [Refs. 5.6]. These procedures were selected to allow 


the simulation to be performed on a personal computer. 


C. TESTING THE BEARING EXTRAPOLATION PROCEDURE UNDER IDEAL 
CONDITIONS | 

Here bearings to the Target were considered without error. Also the bearing rates 
immediately before the turn (used to extrapolate bearings forward in time) and imme- 
diately after the turn (used to extrapolate bearings backward in time) were computed 
exactly with no errors introduced. The only source of error in this case was the linear 
bearing extrapolation. Thus case 1s intended to be an optimistic bound on the accuracy 
of the bearing extrapolation procedure when linear regression 1s used. 

The measure of effectiveness (MOE) for the accuracy of this procedure was the 
percentage of times that the estimated error was within 10% for range, 20° for course, 


and 10°. for speed from the actual parameters of the Target track. These results are 


presented in Table | and Table 2 for lead-lag and lag-lead maneuvers respectively. 
Each number in these Tables resulted from the examination of 18 simulation exper- 
iments (six Target courses and three Target speeds). For example, when the initial dis- 
tance was large, lead and lag angles 70°, and Own Ship speed 15 knts, then five of the 
18 simulations experiments (or 27.8%) resulted in the calculated Target speed being 
within 10% of the actual Target speed. Note that without bearing error, the simulation 


experiments are deterministic. 


Table 1. ACCURACY IN IDEAL CONDITIONS FOR BEARINGS EXTRAPO- 
LATION PROCEDURE (LEAD-LAG MANEUVER) 
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Table 2. ACCURACY IN IDEAL CONDITIONS FOR BEARINGS EXTRAPO- 
LATION PROCEDURE (LAG-LEAD MANEUVER) 


Percentage of times the absolute error 1s 
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As can be seen in Table I and Table 2, the following results were obtained under 


ideal conditions: 


1. Range estimation is good, but only for medium or large initial distances. 


ty 


Accuracy improves with increases in Own Ship speed. 
Course and speed estimation 1s less accurate than range estimation. 


Accuracy improves with increases in lead and lag angle. 


in LS we 


Lead-lag maneuvers give better range estimates than do lag-lead, but not neces- 
sarily better course and speed estimates. 
It is felt that the procedure performed poorly for short initial distances, in part be- 


cause the bearing versus time curve 1s more nonlinear in this case. 


D. TESTING WITH BEARING ERRORS 

For these experiments, the program was run with bearing errors added to the simu- 
lated bearings. Also the bearing rate for each leg was computed using a linear least 
squares technique. 

Only lead or lag angles of 90° were considered, since it was determined in the pre- 
vious step that this is the best course for Own Ship to follow. Accuracy in course and 
speed were not considered since the errors for these parameters were too large, even in 
ideal conditions. An Own Ship speed of 5 knots was not considered for the same reason. 

The MOE for these experiments was the percentage of times (runs) that the esti- 
mated range at time ¢,, was within 5%, 10%. and 20% of the actual range. Each number 
in Table 3 and Table 4 results from 1,800 trials (six Target courses, three Target speeds 


and 100 replications per (course. speed) combination). 


Table 3. ACCURACY IN REAL CONDITIONS FOR BEARING EXTRAPO- 
LATION PROCEDURE (LEAD-LAG MANEUVER) 
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Table 4. ACCURACY IN REAL CONDITIONS FOR BEARING EXTRAPO- 
LATION PROCEDURE (LAG-LEAD MANEUVER) 


Percentage of times the absolute error 1s 
Smaller than 5%, 10%, and 20% of ac- 
tual for range. 
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From Table 3 and Table 4, the following results are obtained: 
]. The higher the Own Ship speed, the better the accuracy. 
2. The best accuracy is obtained for medium initial distances. 
3. Ihe smaller the bearing error; the better thesaceunicy 


4. A lead-lag maneuver gives better accuracy than lag-lead maneuver. 


It is noted that orthogonal regression was also tried on the estimated Target posi- 
tions, but with results simular to those obtained with linear regression. Subsequent ex- 


periments used onlv linear regression. 


E. ANALYSIS OF THE RESULTS 


It appears that even in the absence of bearing errors, the Bearing Extrapolation 
Procedure gives poor estimates of Target course and speed. And when bearing errors 
are introduced, the range errors also become unacceptable. It is the purpose of this 
section to determine why these poor results were obtained. 

1. Differences Between Lead-Leg and Lag-Leg-Estimation 

Figures 2 to 5S show actual and estimated Target positions for four encounter 
geometries with no error in bearing. So the only errors introduced are due to the linear 
extrapolation of bearings. It 1s clear that lag-leg estimation is much more accurate than 
the lead-leg-estimation. It is also clear that any attempt to regress estimated Target 
positions over both legs can result in large errors. 
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Figure 2. Actual (.) and Estimated (*) Target Track for Lead-Lag Maneuver. 
(Target Course = 060°) 
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Figure 3. Actual (.) and Estimated (*) Target Track for Lead-Lag Maneuver. 
(Target Course = 150°) 
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Figure 4. Actual (.) and Estimated (*) Target Track for Lag-Lead Maneuver. 
(Target Course = 030°) 
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Figure 5. Actual (.) and Estimated (*) Target Track for Lag-Lead Maneuver. 
(Target Course = 120°) 


2. Explanation of Why Lag-Leg Estimation is Better than Lead-Leg Estimation 

To help explain why the lag-leg estimation was better than lead-leg estimation, 
scatter plots of bearings versus time for numerous encounter geometries Were obtained. 
Figures 6, 7, and 8 are typical. In Figure 6, a lead-lag encounter is represented. It was 
observed that during the lead-leg (first leg), the bearing rate was generally smaller than 
that of the lag-leg (second leg). Also, and perhaps more importantly, the change in 
bearing rate (bearing acceleration) was also smaller in the lead-leg. When this happens, 
it is possible to use linear regression to extrapolate the lead-leg bearings accurately into 
the future. This makes the lag-leg estimation more accurate. It is clear from Figure 6 
that if linear regression were used to extrapolate the lag-leg bearings into the past, sig- 
nificant errors would result, and these errors would likely cause the lead-leg estimation 
Fovee tess accurate. 

In the same fashion, Figure 7 shows a plot of bearings versus time for a typical 
lag-lead encounter. Here also. the lead-leg (second leg) has a smaller bearing rate and 
bearing acceleration. So linear regression can be used to extrapolate accurately the 


lead-leg bearings backward in time, again making the lag-leg estimation most accurate. 


Unfortunately, however, there are encounter geometries where neither the lead 
nor the lag-legs show small bearing acceleration. For example, Figure 8 shows a case 
where the Target is at high speed, short initial distance, and on course 160°. Here both 
lead and lag-leg bearings are nonlinear, making the Bearing Extrapolation Procedure 
unusable. Extensive experimentation indicates that when all the following conditions 
hold: 

Target course > 150 °, 
Target speed > 5 knots, 


Target range = 40,000 yds; 


then the lead-leg bearings (and the lag-leg bearings) are not significantly linear to be 
accurately extrapolated. and as a result, lag-leg estimation (as well as lead-leg esti- 
mation) 1S poor. 

Finally, it is noted that for very long initial distances (> 40,000 vds), both the 
lead and lag-legs have approximately a linear change in the bearings, making both lead 


and lag-leg estimations reasonably accurate. 
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Figure 6. Bearings in Both Legs for Lead-Lag Maneuver (Target course = 060°) 
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Figure 7. Bearings in Both Legs for Lag-Lead Maneuver (Target Course = 140°) 
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Fieure 8. Bearings in Both Legs for Lead-lag Maneuver at Short Distance and 


Closing Target. 
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F. CONCLUSIONS 


1. 


Even when bearings were without error, the Bearing Extrapolation Procedure was 
unsuccessful in estimating Target course and speed. However, range estimates were 
reasonable. 


When bearing errors were introduced, range estimates also became inaccurate. 


For most encounter geometries, lag-leg estimation gave more accurate results than 
did lead-leg estimation. However. for some important cases both lag and lead-leg 
estimations performed poorly. 


i 


II. POSSIBLE IMPROVEMENTS, AND ANALYSIS OF THEIR 
ACCURACIES 


The lack of accuracy of the proposed procedure, even under ideal conditions, led to 
attempts to improve its performance. The following are some suggested minor changes 
and two suggested major changes that were found to improve the accuracy of the Bear- 


ing Extrapolation Procedure. 


A. MINOR CHANGES 
1. Increase in Leg Length 
In an attempt to improve performance, several different leg lengths were tried. 
It was discovered that leg lengths of 8 to 12 minutes gave the best results. Shorter legs 
gave too little data and longer legs resulted in nonlinear changes in lead-leg bearings. 
Based on these tests, ten minutes was selected as the leg length for subsequent simu- 
lation experiments. 
2. Use of a Quadratic Model to Fair the Lag-Leg 
Another important change from the original simulation that was tried was to 
use a quadratic model to fair the bearings in the lag-leg. This produced better accuracy 
in the estimated Target track parameters. It is noted that current manual TMA methods 
(in particular, the Time-Bearing Plot) use linear fairing. 
3. Increase in Turn Time 
Contrary to Ekelund ranging, if an instantaneous turn 1s considered, large errors 
in the estimated Target track parameters are generated. This 1s because some actual and 
corresponding extrapolated bearings are too close to each other (1. e., a small baseline). 
The intersection of those bearings often results in a estimated Target position far away 
from the true position of the Target. 
Based on these considerations, it was determined that the time required for Own 
Ship to turn should be increased. After many trials, 11 Was found that the best turn du- 
ration is ten minutes, if leg length is also ten minutes. For longer times, the bearing rate 
in the lead-leg does not remain low and constant, reducing the accuracy of the lead-leg 


bearing extrapolation. 


B. ELIMINATION OF THE LEAD-LEG-ESTIMATION 

As shown in Chapter IJ, when linear extrapolation of the bearings from the lag-leg 
is used, increased errors are generated in the estimation of the Target track parameters. 
This occurs because the lag-leg generates a high and variable bearing rate compared to 
the lead-leg. A simple solution to the problem 1s to drop the estimated Target positions 
generated by the extrapolation of the lag-leg bearings. 

One problem with this idea is that a lead-lag maneuver will then produce more ac- 
curate estimates of Target position at the end of the maneuver than will a lag-lead ma- 
neuver. This occurs since in the lead-lag case, the estimated Target positions are 
obtained for times during the second (1.e., the lag) leg. And in the lag-lead case, the es- 
timated Target positions are for times in the first (again, the lag) leg. So to use a lag-lead 
maneuver to estimate Target positions at the end of the maneuver requires that these 
positions be extrapolated in time using derived estimates of Target course and speed. 
This leads to accumulated errors in the final Target position. 

To test for the possible improvements that might result when using only lag-leg es- 
timation. simulation experiments Were conducted with the following parameters: 


1. Speed of Own Ship: 5, 10) andweaawers 


2. Speed cf the Target: 3.) loeand22 seo. 

3. Initia! range: large (= 60.000), medium (= 30,000), and short (= 10,000) vds. 
4d. Only the lead-lag maneuver is considered. 

5. Only the initial lead angle of 90° 1s cosidered. 


G. Bee lenethn: 74 minutes 

7. Vamesserveen legs: 6 minutes: 

§. Standard deviation bearing error 0) ss sal) ace le 

9. Bearings for the first leg were faired and extrapolated linearly. 


10. Bearings in the second leg were faired linearly. 


It is important to note that other values for lead angle were not considered because, 
as discussed in Chapter I], it was found that the best results were obtained for a 90° lead 
angle. 

In order to show how the minor changes improve the performance of this major 
suggested change (i. ¢., elimination of the lead-leg estimation), the simulation was run 


for three different cases: 


1. Without Minor Changes 

The simulation program was run 100 times for each combination of the above 
parameter values. Jt was determined that for the ideal case (0° bearing errors) the ac- 
curacy of the predicted Target parameters was adequate, but with bearing errors intro- 
duced, large errors in all three Target parameters were produced. So it was concluded 
that simply eliminating the lead-leg estimation was not a sufficient improvement. 

2. With Increased Leg Length and Lag-Leg Faired with a Quadratic Model 

Here leg length was increased to ten minutes and the lag-leg was faired with a 
quadratic model. The same simulation experiments were run as before. Only the results 
for lead-lag maneuvers are presented in Table 5. Lag-lead maneuvers gave uniformly 
worse results for reasons discussed above. Additionally, and more surprisingly, the ac- 
curacy for the lag-lead maneuver was not as good as the lead-lag, even when estimating 
the Target parameters at the end of the first leg. This was a confirmation that, in gen- 
eral, the lead-lag maneuver gives better results than do lag-lead [Ref. 7: pp. 3]. 

Results for estimated Target course and speed were very precise only when the 
bearing error was zero and when Own Ship speed was greater than ten knots. When 
errors in bearings Were included, the accuracy in estimated Target course and speed was 
tremendously reduced. For example, for large initial distance and 15 knots Own Ship 
speed the following was obtained: 


1. With no bearing error, the error in estimated speed was alwavs less than 10% and 
the error in course was less than 10° in 83% of the different geometries. 


to 


With a small bearing error (standard deviation = 0.5°), onlv 33.4% of the runs had 
an error in speed of less than 10% . In 15.1% of the runs. the error in course was 
less than 10°. 

One of the conclusions of this study is that the Bearing Extrapolation Procedure, even 
as modified here, does not produce accurate estimates of Target course and speed. With 
Target range. however, more success was obtained. Table 5 shows the percentage of 
simulation runs that produced an estimated Target range at the end of the maneuver 


Within 5°0. 10%, and 20% of the true range. 


Table 5. ACCURACY IN ESTIMATED RANGE WHEN ONLY 
LAG-LEG-ESTIMATION IS USED. (LEG LENGTH = 10 MIN). 


; Percentage of times the absolute error 1s 
Own Ship St. Dev. Initial dis- | smaller than 5%, 10%, and 20% of ac- 
speed bearing er- tance tual for range. 


(Knts) ror 
20 





From Table 5 it can be seen that the larger the Own Ship speed , the betteras 
the estimate of range. This is shown in Figure 9 which is the plot of the accuracy in 


estimated range versus standard deviation of bearing errors when the Own Ship speed 


is 10 and 15 knots. This accuracy is the percentage of times that the absolute error is 


within 20° of the actual for range at the end of the maneuver. 


OWN SHIP SPEED @ 15 Knits. OWN SHIP SPEED = 10 Knts. 





Figure 9. Comparison of Accuracy for Different Own Ship Speeds by Using only 
the Lag-Leg-Estimation. 


Figure 10 shows the comparison of the results of these suggested improvements with the 
results obtained in Chapter II for the central idea. It appears that for short initial dis- 
tances the accuracy of the procedure decreases, but for large and medium distances, even 
with large bearing error (standard deviation = 1.5”), the accuracy increases consider- 


ably. Note that the accuracy for short distances was not adequate in either case. 


Ze 


OWN SHIP SPEED @ 15 Knts. FIRST IMPROVEMENT OWN SHIP SPEED © 15 Knts. USING CENTRAL IDEA 





Figure 10. | Comparison of the Accuracy Between Results by Using only the Lag- 
Leg-Estimation and the Results by Using the Central Idea. 


3. With Previous Minor Changes Plus an Increase in Turn Duration 
In addition to the two minor changes made before, here the turn duration was 
increased to ten minutes and the simulation experiments were repeated. Table 6 shows 


the results for this case. 
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Table 6. ACCURACY IN ESIPIVEATED RANGE WHEN ONLY 
LAG-LEG-ESTIMATION IS USED. (LEG LENGTH = 10 MIN. 
TURN TIME = 10 MIN.) 
















Percentage of times the absolute error 1s 
smaller than %, 10%, and 20% of ac- 
tual for range. 


medi 
shor 
: 
05° 
[target 
Ls 
ea A 

10 
as’ 
ee 







Stee 
bearing er- 
ror 






Own Ship 
speed 
(knts) 






















Initial dis- 
tance 

















Gs 
re 


or 
——> 
i ore) 


2 | to 
‘ 


u 
RS 


By comparing the results obtained for the accuracy in range in Table 5 (Own 
ship turn = 6 minutes) and Table 6 (Own ship turn = 10 minutes), It 1s appears that 
small improvements in precision are obtained by using ten minutes for the turn. 


Figure I] is a graphic representation of this fact. 
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Figure 11. | Comparison of Accuracy for Different Own Ship Turn Times by Using 
only the Lap-Leg-Estimation. 


a. Effect of Target Courses Greater than 150° when Lag-Leg Estimation and 

Minor Changes are Used 

It was stated in Chapter I] that the bearing rate for the lead-legs does not 
remain constant when the Target course is greater than 150° (initial angle on the bow 
less than 30°). This is because range during the maneuver decreases to the point where 
the lead-leg bearings do not maintain linearity. To determine the effect on the esti- 
mation of Target range that this behavior produces, simulations were conducted as be- 
fore except that results were obtained for 12 separate Target courses between 030° and 
175°. Previous results were averaged over Target course, so the effect of a particular 
Target course was not observed. 

The results of one simulation experiment are depicted in Figure 12 and 
presented in Table 7. The parameters used to obtain these results were the following: 


1. Initial distance: 30,000 vds. 


2. Target couses: 030°, 050°, 070°, 090°, 110°, 1307, 150°, 155, 160°, 165°, 170 sae 
ies 


3. larget speeds, 2s lo. and 5 knots. 
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4. Own Ship speed: 15 knots. 
5. Standard deviation of bearing error: 1.5°. 


6. Leg length: 10 minutes. 


7. Turn time: 10 minutes. 





Figure 12. Fraction of Times Absolute Error is within 20% of Actual for Range, for 


Different Target Courses and Speeds. 


a 


Table 7. FRACTION OF TIMES ABSOLUTE ERROR IS WITHIN 20% OF 
ACTUAL FOR RANGE, FOR DIFFERENT TARGET COURSES AND 
SPEEDS: 











030 | 050 | 070] 090 | 110] 130 150 155] 160] 165 | 170 | 175 


15 ose} 0.95] 0.95] 0.97f 0.99} 1.00) 0.94) 0.62] 0.38] 0.26] 0.13] 0.08 
Ps Loaf o.sef o.9af 0.97} 0.95} 0.93) 0.93) o.9af oil 097] oi 095 


From Figure 12 and Table 7. it can be seen that for a Target speed of 5 
knots, approximately equal accuracy was obtained for all Target courses. But for Target 
speeds of 15 and 25 knots, much less accuracy was obtained when the Target course was 
between 150° and 175°. For Target courses less than 150°, there was less observed 


change in the estimation of Target range. 


C. USE OF AN ARCTANGENT MODEL TO EXTRAPOLATE THE LAG-LEG 
BEARINGS 

As presented so far, a major problem with the Bearing Extrapolation procedure 1s 
that the lag-leg bearings cannot be linearly extrapolated forward or backward in time 
while maintaining accuracy. The solution to this problem suggested in the previous se- 
veral sections was simply not to use extrapolated bearings from the lag-leg. The problem 
with this approach is that tracking information is lost when these extrapolated bearings 
are not used. In this section, an attempt will be made to recover these extrapolated 
bearings by using an arctangent, rather than linear, extrapolation model. 

It is known that without change in Own Ship course, bearing to a constant course 
and speed moving Target 1s of the form 


(4 
B, = By + tan™'[ = (t ~ t)], (3.0) 
cpa 





where 


B, 1s the bearing to the Target at CPA2, 


2 Closest Point of Approach 


/, 1s the relative speed between Target and Own Ship, 
Ri q 1S the range to the Target at CPA, 


eis the time at which CPA is to occur. [Ref. 8: pp. 2] 


Assuming that the extrapolated bearings for the lag-leg follow this arctangent 
model, two different attempts were made in order to include these bearings in the sol- 
ution of the problem. 

1. Fitting the Arctangent Model by Using Least Squares 
; An attempt was made to use a least squares technique to find the values of By, 
——. and 1, in (3.0) which produced the best fit to the observed bearing data. Although 
this approach Was ultimately unsuccessful, there may be another least squares procedure 
that will prove more robust. To help any subsequent researcher, the details of this at- 
tempt are presented to illustrate what has already been examined. 
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Now make the change of variable 
}, = tan(B, — Bo), 
or 
Y= a(t — bo). 


Then applying least squares, it is necssary to minimize 
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A necessary condition 1s 
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Dividing (3.3) by 7 and solving for a, 
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Because the determination of Y, for all r, requires knowledge of the value of the 
bearing at CPA (8,), an iterative procedure was attempted where the first value for B, 


was guessed, then with equations (3.5) and (3.6), « and z, were computed. Then by using 


= oes Yn —tan™ [a(t — %)]}, 


{=} 


a new value for B, was obtained. This iteration was continued in an attempt to converge 
on the true B,. Unfortunately for even small errors in bearings, (standard deviation = 
0.5°), the procedure would not reliably converge. 
2. Estimation of the Parameters for the Arctangent Model by Solving a System of 
Linear Equations 
It is known that with three different bearings to a Target that is moving with a 


constant course and speed, the unknown parameters of the equation 
B,= By + tan” [a(t — t)] 


can be determined. [Ref. 7: pp. 2] 
wow let B,, B,, and B, be the bearings obtained at times 4,, 7,, and t, respec- 


tively. Then it can be written 
B, = By + tan” [a(t — 4], 
By = By + tan™'[a(t, — t)], 
ees 18 ten Pele, = 1 Ip 
An attempt will be made to solve for B. «, and 4. Substracting, 
ey Jb tan [a(t = 1g) le tan” [a(t = igle 
B, — By =tan™ '[a(t, — &)] — tan” '[a(t, — &)]). 
PB tanenloi 7h tannic tc 


Now taking the tangent of both sides and knowing the trigonometric identity for the 


tangent of the difference of two angles, the following can be obtained: 


a(t, — 0) 


nr 
L + a°(ly = Io )(t) = fo) 


ro z a(f, — f)) 
Lane ce! > | a, ne ae 
1+ a (03 — fo)t — 4) 


tan(B, — Bs) o:(f3 — fy) 
an eS) ee 


From these follows that 


(2) — fo)(4, — &) tan(B, — B Jage-y- Se (3.7) 
B,-—B 
(5 = tt = 4) tan, = By) = (y= 4) - SP, (3.8) 
l tan(B; — By) 
(13 — fy )(fa — fo) tan(B; — By) == (3 -— 4) - ee eee (3.9) 


Dividing (3.7) by (3.8), (3.7) by (3.9), and (3.8) by (3.9), allows the solution for 4, a, and 
Bee 


| 2 tan(B; ~ By) tan(B = Bile = BG = 4) = (5 = Wie = AN] + (6 = bY 
; (7) — 2) tan(B, — B,)[(3 — 4) tan(B; — By) — (43 — &) tan(B; — B,)] 
tan( B, — B,) tan(B; — By) + B{y - ie tan(B; — B,) tan(B; — B,) 
— (2 — 4) tan(B; — B,)[(3 — 4) tan(B, — By) — (y — 4) tan(B; — B,)] 


(7, — f,) tan(B, — B,) tan(B,; — B,) 


Eo TAME, BDF y= Now OMS 


By = B, - tan Loz, = io)a 


A subroutine to compute these values was written and the simulation of the 
problem was done with the following conditions: 


|. Speed of Own Ship: S37 lOVaneesaxins 


2. Speed of the Tareet> 5, 1a ange 2eenmrs 

3. Initial range: 10,000, 30,000, and 60,000 yds. 
4. Initial angle in the bow for Own Ship: 90°. 
a beeselemarh mC mites: 


6. Time Detwicen jecs; 6 minuies 

7.. Standard deviation bearme emer 0 902 7.0 | andes 

8. Onlv lead-lag maneuver 1s considered. 

9. Bearings for the first leg were faired and extrapolated linearly. 

10. Bearings in the second leg were faired by using a quadratic model. 


11. The extrapolation of the bearings for the second leg was done by using the for- 
mulas obtained above and by using faired bearings from that leg. This was done 


ag 


because the use of raw bearings produced too much error in the estimation of the 
parameters of the Target. 


Results are presented in Table 8. It is emphasized that these results are only 
for lead-lag maneuvers. 


Table 8. ACCURACY IN ESTIMATED RANGE BY USING’ THE 
ARCTANGENT EXTRAPOLATION MODEL FOR THE LAG-LEG 
(LEAD-LAG MANEUVER). 
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Comparing these results with those on Table 3 on page 13, it is clear that this 
procedure performs better than the original (unimproved) Bearing Extrapolation Proce- 
dure. However, the results are about the same as those in Table 8, indicating that sim- 
ply not using the lag-leg extrapolated bearings is just about as effective as trying to 


extrapolate them with the arctangent model. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


fe CONCLUSIONS 


. The proposed procedure does not allow estimation, with adequate accuracy, of the 


Target track parameters when linear bearing extrapolation is used for both legs 
(central idea). Linear bearing extrapolation is only valid for lead-legs where low 
and constant bearing rate 1s obtained. For lag-legs the bearing rate is generally 
higher and more variable. 


Adequate accuracy for course and speed is not obtained using any of the examined 
procedures unless bearing error is zero (ideal conditions). 


Adequate accuracy in estimated range for large and medium initial distances can 
be obtained 1f: 


a. Quadratic faired bearings from the lag-leg are used in combination with linear 
extrapolated bearings from the lead leg. 


T 


Lead-lag maneuver 1s performed. 


Lead angle close to 90° 1s used. 


a 2 


Leg length 1s larger than seven minutes and less than twelve minutes. 


& 


Own ship turn time ts larger than six minutes and less than ten minutes. 


f. Angle on the bow is not smaller than 30°. If it is, then Target speed is not 
greater than 5 knots. 


Adequate accuracy in range for short initial distances is obtained only when the 
Target is not closing the Own Ship. 


2. RECOMMENDATIONS 


If this procedure is implemented, use only extrapolated bearings from the lead leg 
(low and constant bearing rate) and faired bearings for the lag-leg (high and vari- 
able bearing rate). 


Further research should be done to obtain a model that allows use of the extrapo- 
lated bearings obtained from the lag-leg. 


APPENDIX A. DISCUSSION OF SELECTED SOLUTION METHODS 


A. SPIESS RANGING 

Spiess theory establishes that the problem of obtaining satisfactory solutions for a 
bearings-onlv approach against targets moving with a constant course and speed, can 
be solved by using four bearings and at least one change of course or speed by Own 
Ship. With this in mind, several graphical methods capable of giving quick reliable sol- 
utions Were developed. [Ref. 3] 

From the Passive Ranging Manual, Volume III, the basic graphical method for 
Spiess ranging 1s explained as follows: 

Given three bearings to the Target, observed at times 4,, #, , and 4, together 
with a fourth time, 7,. the locus, L, of all possible Target positions at time, 4, for all 
possible Target tracks which satisfy the three bearings at 4, 4, and 4%, and which 
maintain a constant course and speed. is a straight line. The actual position of the 


Target at time, ¢,. 1s determined by the intersection of L with the observed bearing 
at time, 4. (provided SSK changes course and speed). 


This locus, L, can be determined by picking out two arbitrary Target tracks 
which satisfy the three-bearing conditions. and plotting the position of the Target 
on these tracks at time 4, Then L is the line through these two plotted positions. 
This is the basis of the Spiess four-bearing TMA. The bearing lines at 4, 1, , and 4, 
are used to find the locus of Target positions at z,. The intersection of this locus 
with the bearing line at 4, determines one point on the track. Another point 1s de- 
ternuned by using the bearing lines at 4. 4, and 4, to find the locus of Target posi- 
tions at 4. The intersection of this locus with the bearing line at 4, determines 
another point on the track. The track is then the hne joining the two constructed 
Target positions. 


Because bearings are not precise enough to apply this method, faired bearings are 


usually used. 


B. EKELUND RANGING 

The basic Ekelund solution for a constant course and speed moving Target is based 
on a two leg maneuver by Own Ship. During these legs, Own Ship records actual] Target 
bearings. The rate of change of the bearings is computed for each leg bv applving linear 
least squares techniques. This method assumes that during the turn the distance that 
Own Ship moves and the change of angle in the bow can be disregarded; 1.e., assumes 
an instantaneous turn. Given this assumption. the range to the Target can be computed. 


Using Figure 13. derivation of range is as follows: 


we 
OO 


¢— target track 


<——— own ship track <———leg two 





Figure 13. Geometry for Ekelund Ranging Method. 


The velocity of the Target across the line of sound for Own Ship’s first leg is 


Op , 
Van = Vy sind, = R, aoe = RB (A.1) 
and for leg 2 
: op . 
Van = Vp sin 8, = Ry ee = R,Bp (A.2) 


Own ship’s velocity across the line of sound for leg 1 can be expressed as 


| co 
Vao1 = Vo, Sind, = R a = Ri Bo (4.3) 
and for leg 2 
| 60, 
Vigor = Voz SIN G2 = Ry Aa oy Ry Bo? (A.4) 


where 
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B,, = Bearing rate of Target in leg i. 
R, = Range to the Target at the turn. 
Kk, 


= Range to the Target after the turn. 


V,, = Velocity of Own Ship in leg 1. 
l= Velocity ofthe Tarcemummles 


Assuming Target velocity remains constant during Own Ship maneuvering; te., 
V,=V,= V, and that the angle in the bow also remains approximately constant; 1.e., 
0,=@,=80, and that the range to the Target does not appreciably changes, sie 
R, = Rk, = R,, then equations (A.1), (A.2), (A.3), and (A.4) can be written as 


SiG =e (A.5) 
I’, sin @ = R,Bp, (A.6) 
V,, sind, = RBar, (A.7) 
V. sin @ = R,B,». (A.8) 


and by subtracting (A.7) from (A.5) and (A.8) from (A.6) 


Vy sin@ = (sing, = te =a eee 


Vesmn@— V5 sin 6, = R{ Bp — Boo) = RB. 
Subtracting these last two equations vields 


74, sin @, — Fy sin @9 


k= 
° B, — B, 


% 


where, B, = bearing rate of Target relatrve to Own Ship in leg 1. 


The problems in accuracy with this method are: 


1. It assumes that Own Ship can perform the two legs maneuver with an instantane- 
OUS TUihINMel. 2 pp. 1-3] 


2. The Ekelund Range Equation is dependent on the ratio of the bearing rates devel- 


oped on the two legs of the rangmg maneuver which are small and not very precise 
VAs) | Clam (onesie): 
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3. In practice, measuring the bearing rate involves a bearing smoothing process which 
takes place over a portion of each leg and results in bearing rates at times signif- 
icantly different from the time of the turn [Ref. 9: pp. 3-2]. 

C. CHURN METHOD 

The basic concept of CHURN TMA is similar to that of the strip plot and includes 
in its solution the Spiess TMA for the particular case where there are only four bearing 
lines computed (Ref. 9: pp. 7-1, 7-3]. Figure 14, depicts the regular CHURN TMA 
method. 


estimated 
target track 


Bi = bearing 
at time | 


<j 


own ship track 





Figure 14. Geometry for CHURN Method. 


Own Ship performs a two leg maneuver during which it measures bearings to a 
Target that 1s moving with constant course and speed. The CHURN method takes those 
bearings and fits a Target track to them using the perpendicular distance between a given 
bearing line and the corresponding Target position on the fitted track as a measure of 
the goodness of fit. Applving the concept of least squares, the best Target track is that 


which then minimizes the sum of squares of these distances.[Ref. 9: pp. 7-2 | 
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APPENDIX B. VARIABLES OF THE SIMULATION PROGRAM 


. A,B.C,D.E.F. Entries of the A matrix for quadratic regression in the first lem 
These values were computed by using APL [Ref. 10]. Then they were introduced 
as initial data. 


AQ.A1,A2. Coefficients for the quadratic model in the first leg. 
AOLAI!!A22. Coefficients for the quadratic modeliin the Second leg, 


B22. Temporary variable used to compute the estimated position of the target on 
the estimated target track at the begining of the maneuver. 


BEERR» Beating enor 
. BEEXTRUN). Vector of extrapolated bearings. 


. BEF. Temporary variable used to compute the estimated position of the target on 
the estimated target track at the end of the maneuver. 


. BECN). Vector of bearings to the target. 

. BEIIN. Y intercept for the linear model of bearings versus time in the first leg. 

. BEI(IN). Temporarv vector of bearings to the target. 

. BE2IN. Y intercept for the linear model of bearings versus time im the seconds 
. BMEAN. Mean of the bearing to the target. Same variable is used for both legs. 
. BO. Bearing to the tareet am en 

. CC. Counter used in the computation of the precision of the procedure. 

. CINCO. Times that absolute value of the range error is less than 5%. 

. CINCO. Times that absolute value of the speedtermorsicslcsoumn aimee c 

. CINCO2. Times that absolute valtievor tie course error 1s less staciamee 

. COFL(2). Vector of own ship courses for the first leg. 

. COST. Own ship course for the’ secondsies 


. C5. Percentage of times that error in estimated target course is Jess than 5°. 





. C10. Percentage of times that error in estimated target courseus less than 103 


. C20, Percentage of times that error in estiniated (anges courses less tiie) 





DEGMIT Seive presseduim degree. 

. DIEZ. Times that absolute value Of the tangeterroris leccune ama 7 

. DIEZI. Times that absolute value of the speed error is less than 10%. 
. DIEZ2. Times that absolute value of the course error is less than 10°. 


. DSEED. Seed for the pseudorandom generation of normal bearing errors. 


28. ECTLR. Estimated target course obtained by linear regression on the target fixes. 
29. ERANTL. Estimated range to the target at the end of the maneuver. 

30. ESLOPE. Temporary variable corresponding to the target course. 

feo tL. Estimated speed of the target. 

32. EXTUN). Vector of estimated target positions in the X axis. 

Semele t. Y intercept of the estimated target course. 

34. EX TUN). Vector of estimated target positions in the Y axis. 

Seb e(]N): Vector of faired target bearings. 

36. IN. Time to complete the maneuver in 1/3 of minutes. 


37. INFOUN). Flag used to signal when Estimated X position of the target is equal 
to the X position of the own ship. 


38. IA(3). Vector of initial target distances from own ship in the X axis. 
39. IN(3). Vector of initial target distances from own ship in the Y axis. 


40. LLL. Counter to determine the amount of fixes used in the regression to determine 
target track. 


41. PER. Percentage error in range computed from the estimated and the simulated 
famee tO the target at the end of the maneuver. 





42. PERI. Percentage error in speed computed from the estimated and the simulated 
Farcet speed. 


43. PER2. Amount of error in course computed from the estimated and the simulated 
target course. 


feeer il. Piin radians. 


45. RS.T.UV.W. Entries of the A-' matrix for quadratic regression in the second leg. 
These values were computed by using APL [Ref. 10]. Then they introduced as in1- 
tral data. 


46. RANGEF. Simulated range to the target at the end of the maneuver. 


47. RCBFLI. Rate of change of bearing during the first leg. This value 1s computed 
assunung the linear model for the bearings Versus time. 


48. RCBSLI. Rate of change of bearing during the Second leg. This value is computed 
assuming the linear model for the bearings versus time. 


49. RS. Percentage of times that error in estimated range to the target 1s less than 5%. 


50. R10. Percentage of times that error in estimated range to the target is less than 


OP 





51. R20. Percentage of times that error in estimated range to the target is less than 





52. SC1(6). Vector of simulated target courses. 


33. SD. Standard deviation of bearing errors. 


54. 


55; 


56. 
ee 


58. 
Bye 
60. 


. SUMB2. Summation of the square of the bearing values. 


68. 
69, 
dO: 


7 


mee 


SLOPECUIN). Vector of faired bearings to the target expressed as slopes in the X-Y 
plane. 


SLOPEX(UIN). Vector of extrapolated bearings to the target expressed as slopes 
in X-Y plane. 


SMALL. Constant that defines equality between simulated Own ship position and 
estimated target position. 


SMALLB. Constant used to determine when the faired or extrapolated target 
bearing correspond to a value of 90° or 270°. 


SO(3). Vector of simulated own ship speeds. 
SST(3). Vector of simulated target speeds. 


SUMB. Summation of the bearing values. 


.SUMT. Summation of time. 

. SUMTB. Summation of the products of bearing and time. 

. SUMT2. Summation of the squares of the time. 

. SUMT2B. Summation of the products of square of time and bearing. 
. SUM A. Summation of the estimated target positions in the X axis. 


. SUMAY. Summation of the product of estimated target positions in the X axis 


and the estimated target position in the Y axis. 
SUMX2. Summation of the square of the estimated target positions in the X axis. 
SUMY. Summation of the estimated target positions in the Y axis. 


SUMY2. Summation of the square of the estimated target positions in the Y axis. 


. SNTCUN). Vector of simulated target positions in the X axis. 
12 
fis} 
74, 
TS, 
76. 
ngs 
78. 
a) 
SO. 
oi 
82. 


SYTUN). Wector of simulated target positions in the Y axis. 





S5. Percentage of times that error in estimated target speed is less than 5% 


; 


N 


10, Percentage of times that error in estimated target speed is less than 10% 


Ca) 


20. Percentage of times that error in estimated target speed is less than 20% 


: ak 
: 


WELZ] Vector of leglensitine 

TEMP. Temporary variable. 

Jee os. .Viean ofthe unic 

TO. Time when CPA occurs. 

TWO. 2 expressed as double precision constant. 
VARERB. Variance of bearing error. 


VEINTE. Times that absolute value of the range error is less than 20%%. 


a4 


S23. 
84. 
oS. 
SO. 
87. 
88. 
89. 
90. 
oN. 
Be. 
98. 
94. 
oD. 


VEINII. Times that absolute value of the speed error is less than 20%. 
VEINT2. Times that absolute value of the course error is less than 20°. 
VD. Ratio between relative speed of target and own ship and range at CPA. 
AO(UIN). Vector of own ship positions in the X axis. 

NOENTRUN). Vector of extrapolated own ship positions in the X axis. 
AMEAN. Mean of target positions (fixes) in the X axis. 

AFF. Estimated target position at the end of the maneuver in the X axis. 


X22. Estimated target position at the begining of the maneuver in the X axis. 





YFF. Estimated target position at the end of the maneuver in the Y axis. 
YMEAN. Mean of target positions (fixes) in the Y axis. 

YOUN). Vector of own ship positions in the Y axis. 

YOEATR(UN). Vector of extrapolated own ship positions. 


Y22. Estimated target position at the beggining of the maneuver in the Y axis. 
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APPENDIX C. SIMULATION PROGRAM CODE 
PROGRAM SIMTMA 


evevevevevelededeteselelotetedededetisesedesedesedesededesedtedtedesesedescKdedesesiesedesedidesesededtesedededtesesesesedeiedededensededtedehkexe 


ale 


THIS PROGRAM SIMULATES A PASSIVE BEARINGS-ONLY TARGET MOTION ANALYSIS* 


PROCEDURE WHICH IS PERFORMED BY A SUBMARINE (OWN SHIP) AGAINST A “ 
TARGET WHICH IS MOVING WITH CONSTANT COURSE AND SPEED. % 
THE OWN SHIP PERFORMS A TWO LEG MANEUVER. EACH OF THEM WITH CONSTANT * 
COURSE AND SPEED. a 
{le fa he et we ue cits ris 

we 

WRITTEN BY CMR. CARRERO CUBEROS, BERNABE * 

VENEZUELAN NAVY ¥ 

rhs 
Sesevevevesevevevesedesevedesesedesesedesese Seve secede Fete sete cece See Tee SoS De Ve SEDO S CLOSET Se eee SCTE ERC TC UC DC TR TR RTO UN ae eae 


PARAMETER(IN = 78) 

INTEGER TEL(2), INFOC(IN), LLL, CC, CINCO, DIEZ, VEINTE, 
&CINCO1, DIEZ1, VEINT1, CINCO2, DIEZ2, VEINT2 

DOUBLE PRECISION IX(3), IY(3), DSEED, COFL(2), COSL, SCT(6), SD, 
&SST(3), VARERB, BEERR, EXT(IN), EYT(IN), SO(3), SXT(IN), SYTCIN), 
&XO(IN), YOCIN), BE(IN), SLOPE(IN), BEEXTR(IN), XOEXTR(IN), 
&YOEXTR(IN), SLOPEX(IN), PI, SMALL, DEGPI, TWO, TEMP, SUMX, SUMX2, 
&SUMY2, SUMY, SUMXY, XMEAN, YMEAN, ESLOPE, SMALLB, EYINT, ECTLR, 
&B22, BFF, X22, Y22, XFF, YFF, ESPETL, ERANTL, RANGEF, SUMT, SUMB, 
&SUMT2, SUMB2, SUMTB, TMEAN, BMEAN, RCBFL1, RCBSL1, BE1(IN), BE1LIN, 
&BE2IN, FBEC(IN), AO, Al, AZ, SUMT2B, A, B, CG, Dy E,eh sy 40u eee 
&A21, R, S, T, U, V, W, BO, TO, VD, PER, PER1, PER2, R5, R10, R20, 
fe5) S10, S20 pmCo cle emalg 

DATA PI/3. 1415926535898D0/, SMALL/0O. 1D0/, DEGPI/180.D0/, TWO/2.D0/, 
&SCT/153. DO, 158. DO, 163. DO, 167. DO, 171,D0. 17510) aint) oo ere 
&$0/30000. DO, 20000. D0, 10000. DO/, SST/50000. DO,30000. DO, 10000. DO/, 
&COFL/088.D0,272.D0/, SMALLB/.00001D0/, A/118.4785317D0/, 
&B/-3.77185762D0/, C/.02947719689D0/, D/. 1205858692D0/, 
&E/-.0009459915779D0/, F/.000007448752582D0/, R/. 3438423645D0/, 
&S/-.0450738916D0/, T/.001231527094D0/, U/.007603190052D0/, 
&V/-. 00023091133D0/, W/.000007448752582D0/, 
&1Y/60000. DO, 30000. D0, 10000. D0/, IX/0.DO,0. DO,0. DO/ 
PRINT**,'ENTER VARIANCE OF BEARING ERROR (SQUARE DEGREES) AND' 
PRINT*:,' INITIAL SEED (7 DIGIT INTEGER) TO BEGIN THE SIMULATION’ 
READ** , VARERB , DSEED 


** TRANSFORMATION OF DATA THAT INVOLVES DEGREES TO RADIANS 
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10 


SD = SQRT( VARERB) 
SD = SD*PI/DEGPI 
DO-s0rle = 1 2 
COFL(I) = COFL(1)*PI/DEGPI 
CONTINUE 
DO; GO 3 16 
SCi( 1 )-=7SCil 1) Ply Deore 
0 CONTINUE 
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SELECTION OF TWO DIFFERENT LENGTH LEG 
MoMZ2C0G l = ine 


SELECTION FROM THREE DIFFERENT INITIAL DISTANCES BETWEEN TARGET AND 
OWN SHIP 


POfI279 MMMM =~ 1.3 


SELECTION FROM THREE DIFFERENT SPEEDS FOR OWN SHIP 
DO 1100 L = 1,3 


* SELECTION FROM LEAD-LAG AND LAG-LEAD MANEUVER AND INITIALIZATION 


DO 950 M = 1,1 
CC = 0 


VEINTE 
CINCO1 
DIEZ1 = 0 

VEINT1 = 0 
CINCO2 = 0 
DIEZ2 = 0 

VENT 2 >= 0 


Ge 1 


SELECTION FROM SIX DIFFERENT TARGET COURSES 
DO 900 J = 1,6 


SeoeeresctlION FROM THREE TARGET SPEEDS 


DO 850 K = 1,3 


ieee tiTION OF THE SIMULATION 100 TIMES WITH DIFFERENT SEED 
BO 7257013 1 100 


INITIALIZATION AND SIMULATION OF TARGET POSITION EVERY 20 SECONDS 
DOME Cel 1o= oi 2 ei) 


Se Wen), DO 
Sheu) = 0.00 
XOC II) = 0. DO 
YOCII) = 0. DO 

BEC II) = 0. DO 
Beat) =a. DO 
FBECII) = 0. DO 
SLOPE CLI) = 0. DO 
SLOPEX( II) = 0. DO 
XOEXTR( II) = 0. DO 
YOEXTR( II) = 0. DO 
ERE OR CID =a Osu 
INFOCII) = 0. DO 


REV (DE LECII)) / 1803 p0 
Ser = 1 OMI) + COSCSCiGd SE IWwO)*SSTCK ) “TEMP 
Se ee — YOM eee SINC SEG) pane) S51. Kk) * TEMP 


150 CONTINUE 


ate 
oA 


ute 
ee 


INITIALIZATION OF REGRESSION TO COMPUTE RATE OF CHANGE OF BEARING IN 
Zool LEG. 

Ser? == 030 

SUMT2 = 0. DO 


SUMB2 = 0. DO 
SUMB = 0. DO 
SUMTB = 0. DO 
SUMT2B = 0. DO 

DO 7/50 N = 1 Tener rs 


%* SIMULATION OF POSITION OF OWN SHIP EVERY 20 SECONDS IN THE FIRST LEG: 
TEMP = (DBLE(N))/180. DO 
XO(N) COS(COFL(M) - PI/TWO)*SO(L)*TEMP 
YOCN) SINCPI/TWO - COFL(M))*SO(L)*TEMP 


* SIMULATION OF THE BEARINGS TO THE TARGET IN THE FIRST LEG. 
CALL NORRN(DSEED , BEERR) 
IF(ABS(SYT(N) - YO(N)) .GE. SMALL)THEN 
IF(ABS(XO(N) - SXT(N)) .GE. SMALL)THEN 
IF(YO(N) . LT. SYT(N))THEN 
IF(XO(N) .LT. SXT(N))THEN 
BE(N) = ATAN((SXT(N) - XO(N))/(SYT(N) - YO(N))) 
ELSE 
TEMP = ATAN((SXT(N) - XO(N))/C(SYT(N) - YO(N))) 
BE(N) = TWO*PI + TEMP 
END IF 
ELSE 
IF(XO(N) .LT. SXT(N))THEN 
TEMP = ATAN((SXT(N) - XO(N))/(SYT(N) - YO(N))) 
BE(N) = PI + TEMP 
ELSE 
TEMP = ATAN((SXT(N) - XO(N))/(SYT(N) - YO(N))) 
BE(N) = PI + TEMP 
END IF 
END IF 
BE(N) = BE(N) + BEERR*SD 
IF(BE(N) .GE. TWO*PI)BE(N) = BE(N) - TWO*PI 
ELSE 
IF(SYT(N) .GT. YO(N))THEN 
BE(N) = 0. D0 
ELSE 
BE(N) = PI 
END IF 
INFO(N) = 99 
END IF 
ELSE 
IF(SXT(N) .GT. XO(N))THEN 
BE(N) = PI/TWO 
ELSE 
BE(N) = 3. DO*PI/TWO 
END IF 
END IF 


* COMPUTATION OF COEFFICIENTS FOR QUADRATIC AND LINEAR LEAST SQUARES 
* FITS FOR THE BEARINGS IN THE FIRST LEG. 
SUMT = SUMT + DBLE(N) 
SUMT2 = SUMT2 + (DBLECN) )**2 
IFCBEC(N) .Gi. 3, 00°P1/ We) aN 
BEI(N) = BE(N) - TWO*PI 
ELSE 


4§ 


750 


af 


BE1(N) = BE(N) 
END IF 
SUMB = SUMB + BE1(N) 
SUMB2 = SUMB2 + BE1(N)**2 
SUMTB = SUMTB + (DBLE(N))*BE1(N) 
SUMT2B = SUMT2B + (DBLE(N)**2)*BE1(N) 


CONTINUE 

TMEAN = SUMT/DBLE( TEL(1I)-9) 

BMEAN = SUMB/DBLE(TEL(1I)-9) 

RCBFL1 = (SUMTB - SUMT*BMEAN)/(SUMT2 - SUMT*TMEAN) 
BE1IN = BMEAN - RCBFL1*TMEAN 

ACI =e SUNB + S<SsuUMIB + T<SUMDZB 

All 
A21 


So SUNB 2 UsSUMIB + V*SUMDZB 
TSUMB + V*SUMTB + W*SUMT2B 


* COMPUTATION OF THE OWN SHIP POSITION AT THE FICTITIOUS TURN 
TEMP = VELECTELCI))7 180: 0 

XO(TEL(I)) = COSCCOFL(M) - PI/TWO)*SO(L)*TEMP 
TOC ELC), 


SINC PI/TWO - COFL(M))*SO( L)*TEMP 


* COMPUTATION OF FAIR BEARINGS DURING THE FIRST LEG USING THE RESULTS 


af 


oe 


* FROM LINEAR OR QUADRATIC FIT. THIS DEPEND IF THE MANEUVER IS LEAD-LAG 
* OR LAG-LEAD. 
PO 755 °N = Irene r 9 


IF (M .EQ. 1)THEN 
FBE(N) = RCBFL1*DBLE(N) + BE1IN 

ELSE 
FBE(N) = AO1 + Al11*DBLE(N) + A21*(DBLE(N)**2) 

END IF 

IF(FBE(N) .LT. 0.DO)FBE(N) = TWO*PI + FBE(N) 

IF((ABS(FBE(N) - PI/TWO) .LT. SMALLB) . OR. 
(ABS(FBE(N) - 3.D0*PI/TWO) .LT. SMALLB))THEN 
SLOPE(N) = 0. DO 

ELSE 
SLOPE(N) = 1.DO0/TAN(FBE(N)) 

END IF 


CONTINUE 


* IN CASE OF LAG-LEAD MANEUVER (M=2), CALL SUBROUTINE TANGENT TO COMPUTE 
getHE COEFFICIENTS THAT RELATES THE BEARINGS TO THE ARCTANGENT MODEL. IN 
* CASE OF LEAD-LAG MANEUVER USE THE LINEAR MODEL ALREADY COMPUTED. WITH 
* THIS INFORMATION, COMPUTE THE EXTRAPOLATED BEARINGS FROM THE FIRST LEG 


PEG nO we 2 )CALLeINGENTCPBE IN ope 10. VD) 
DOP7GC °N — TELCI)+10, 2*1E LCD) 


PAC sora lenin 
BEEXTR(N) = RCBFL1*DBLE(N) + BE1IN 
ELSE 
BEEXTR(N) = BO + ATANCVD*(DBLE(N) - TO)) 
eee 
IF(BEEXTR(N) .LT. 0. DO)BEEXTR(N)=TWO*PI+BEEXTR(N) 
eeCABS( BEEXERON) - Pl7awGjeea!. SMALLB) . OR. 
(ABS(BEEXTR(N) - 3.D0O*PI/TWO) . LT. SMALLB))THEN 
SLOPEX(N) = 0. DO 
ELSE 
SLOPEX(N) = 1. DO/TANCBEEXTR(N) ) 


49 


BND TE 
760 CONTINUE 


* SIMULATION OF EXTRAPOLATED POSITIONS OF OWN SHIP OUT OF THE FIRST LEG. 
DO 765 N = Thing 10, 2 7ELCD) 
TEMP = (DBLE(N))/180. DO 
XOEXTR(N) = COS(COFL(M) - PI/TWO)*SO(L)*TEMP 
YOEXTR(N) = SINCPI/TWO - COFL(M))*SOCL)*TEMP 


765 CONTINUE 


* SIMULATION OF FAIR BEARING TO THE TARGET AT THE TURN, BY USING THE 
* LINEAR MODEL. 
FBECTEL(I)) = RCBFL1*DBLECTEL(I)) + BE1IN 
IF( FBECTEL(1)). ahi. ODO) EBB GTEL( 1) )e= WC ee 
& FRECIELGL)) 
IFC FBECTEL(I)) .GE.” TWO*PI)FBECTEL(1)) = PRECIELG a 
& TWO*PI 


* DETERMINATION OF OWN SHIP COURSE FOR THE SECOND LEG ACCORDING WITH THE 
* LINE OF SIGHT TO°WHE TARGET Alii nent 
ie Chie BOS) tHeN 
COSL = FBEC(TEL(I)) - PI*(.48888889D0) 
ELSE 
COSL = FBEC(TEL(I)) + PI*(. 48888889D0) 
END a 
IFCCOSL .LT. 0.DO)COSL = TwO*PI + COSL 
IF(COSL .GT. TWO*PI)COSL = COSL - TwO*PI 


* INITIALIZATION OF REGRESSION TO COMPUTE RATE OF CHANGE OF BEARING IN 
SECOND LEG. 


sf. 
“nn 


SUMT = 0. DO 
SUMT2 = 0. DO 
SUMb2 = 0. DO 
SUMB = 0. DO 
SUMTB = QO. DO 
SUMT2B = 0. DO 
DO°778 N = TELCI) +106 2.7 Era, 
* SIMULATION OF POSITION OF OWN SHIP EACH 20 SECONDS IN THE SECOND LEG. 
TEMP = (DELECN) = DELECGIELGI oy ceed 
XOCN) XOCTELC1}) #8S0(CL) “TEMP. COSCCOSh=P1/TiGD 
YOCN) YOCTEL(I)) - SOCL)*TEMP*SINCCOSL-PI/TWO) 


3. 


SIMULATION OF THE BEARINGS TO THE TARGET IN THE SECOND LEG. 
CALL NORRN(DSEED, BEERR) 
IFCABSCSYICN) = YOCN)) 2 GE. star THEN 
IFCABSC(CSXTCN) - XOCN)) .GE. SMALL)THEN 
IF( YOUN) . LE. SYICN) tier 
IFCXOCN) . LI. SXT(CN))THEN 
BE(N) = ATANC(SXT(N) - XOCN))/CSYTCN) - YOCN))) 
ELSE 
TEMP = ATANCCSXTCN) - XOCN))/CSYTCN) - YOCN))) 
BECN) = TWO*PI + TEMP 
END IF 
ELSE 


tf" 
o> 
— 


IF(XO(N) .LT. SXT(N))THEN 
TEMP = ATAN((SXT(N) - XO(N))/(SYT(N) - YO(N))) 
BE(N) = PI + TEMP 
ELSE 
TEMP = ATAN((SXT(N) - XO(N))/(SYTCN) - YOCN))) 
BE(N) = PI + TEMP 
END IF 
END IF 
BE(N) = BE(N) + BEERR*SD 
IF(BE(N) .LT. 0.DO0)BE(N) = TWO*PI + BE(N) 
IF(BE(N) .GE. TWO*PI)BE(N) = BE(N) - TWO*PI 
ELSE 
IF(SYT(N) .GT. YO(N))THEN 
BE(N) = 0.DO 
ELSE 
BE(N) = PI 
END IF 
INFO(N) = 99 
END IF 
ELSE 
IF(SXT(N) .GT. XO(N))THEN 
BE(N) = PI/TWO 
ELSE 
BE(N) = 3.DO*PI/TWO 
END IF 
END IF 


* COMPUTATION OF COEFICIENTS FOR LINEAR AND QUADRATIC FITS FOR THE 
* BEARINGS IN THE SECOND LEG 
SUMT = SUMT + DBLE(N) 
SUMT2 = SUMT2 + (DBLE(N) )**2 
PEC EO. THEN 
PEChE CN) eeals ~ ©>DO*E I) EnEN 
BE1(N) = BEC(N) - TWO*PI 


ELSE 
BE1(N) = BE(N) 
EN ue 
ELSE 
BEwChoe— BECN) 
END IF 


SUMB = SUMB + BEI(N) 
SUMB2 = SUMB2 + BE1(N)**2 
SUMTB = SUMTB + (DBLE(N))*BE1(N) 
SUMT2B = SUMT2B + (DBLE(N)**2)*BE1(N) 
778 CONTINUE 
TMEAN = SUMT/DBLE(TEL(1)-9) 
BMEAN = SUMB/DBLE(TEL(1)-9) 
RCBSL1 = (SUMTB - SUMT*BMEAN)/(SUMT2 - SUMT*TMEAN) 
BE2IN = BMEAN - RCBSL1*TMEAN 


AO = A*SUMB + B*SUMTB + C*SUMT2B 
Al = B*SUMB + D*SUMTB + E*SUMT2B 
A2 = C*SUMB + E*SUMTB + F*SUMT2B 


* COMPUTATION OF FAIR BEARINGS DURING THE SECOND LEG USING LINEAR 
sek QUADRATIC FIT. THIS DEPEND ON THE TYPE OF MANEUVER. 


DO 779 N = TEL(1)+10,2*TEL(I) 
IF(M .EQ. 1)THEN 
FBE(N) = AO + Al*DBLE(N) + A2%*(DBLE(N)**2) 
ELSE 
FBE(N) = RCBSL1*DBLE(N) + BE2IN 
END IF 
IF(FBE(N) .LT. 0.DO)FBE(N) = TWO*PI + FBE(N) 
IF((ABS(FBE(N) - PI/TWO) .LT. SMALLB) .OR. 


& (ABS(FBE(N) - 3.D0O*PI/TWO) .LT. SMALLB))THEN 
SLOPE(N) = 0. DO 
ELSE 
SLOPE(N) = 1.D0/TAN(FBE(N)) 
EN Derr 
119 CONTINUE 


* IN CASE OF LEAD-LAG MANEUVER (M=1), CALL SUBROUTINE TANGENT TO COMPUTE 
* THE COEFFICIENTS THAT RELATES THE BEARINGS TO THE ARCTANGENT MODEL. IN 
* CASE OF LAG-LEAD MANEUVER USE THE LINEAR MODEL ALREADY COMPUTED. WITH 
* THIS INFORMATION, COMPUTE THE EXTRAPOLATED BEARINGS FROM THE SECOND LEG 

IF(M . EQ) CALL INGENTCEBEIN SM, BORO yy 1) 

DO 781 N = IQPEEG) 2 

1IFCM . EO] an 
BEEXTR(N) = BO + ATAN(CVD*(DBLE(N) - TO)) 


ELSE 
BEEXTR(N) = RCBSL1*DBLE(N) + BE2IN 
END IF 
# WRITE( 1, '(F9.6,1X,13)' )BEEXTR(N),N 


IF(BEEXTR(N) .LT. 0.DO)BEEXTRC(N)=TWO*PI+BEEXTR(N) 
IF( (ABS( BEEXTR(N) = Pl/ TWO) eal. SMAECE Re OR: 


& (ABS(CBEEXTR(N) - 3.DO0*PI/TWO) . LT. SMALLB) )THEN 
SLOPEX(N) = 0. DO 
ELSE 
SLOPEX(N) = 1. DO/TAN(BEEXTR(N) ) 
END IF 
781 CONTINUE 


* SIMULATION OF EXTRAPOLATED OWN SHIP POSITIONS OUT OF THE SECOND LEG. 
DO-782°N = 1, TEE 1-9 
TEMP = (DBLE(TEL(1))> = DBLEGN oy 130. DC 
XOEXTR(N) = XOCTEE(])) —= COs(GGsk.- Pl, iva) 


& SO( L)*TEMP 
YOEXTR(N) = YOCTEL(I)) + SINCCOSL - PI/TwO)* 
& SO(L)*TEMP 
PO 2 CONTINUE 


* INITIALIZATION OF VARIABLES USED IN THE LINEAR REGRESSION TO COMPUTE 
* COURSE OF THE TARGET. 


SUMX = 0. DO 
SUMK2 = 0. DO 
SUMY2 = 0. DO 
SUMY = 0. DO 
SUMAY = Q. DO 
LLL = 0 


* COMPUTATION OF EXPECTED POSTTION (X,Y) OR THE= faker: 


DOmveeeN = 1. ote ECT) 


* ELIMINATION OF 18 CENTRAL ESTIMATED POSITIONS OF THE TARGET TO ALLOW 
* THE MANEUVER OF THE OWN SHIP IN ORDER TO CHANGE COURSE FROM FIRST LEG 
gato sECOND LEG (SIX MINUTES) 


PhCCONS GT. TEL 1 )-9) . AND: ((NeeIEEe CTELC1)+9))) 


& GO TO 783 
IFCINFO(N) .EQ. 99)THEN 
EXT(N) = XOCN) 
EYT(N) = SLOPEX(N)*(XOC(N) - XOEXTRCN)) + YOEXTR(N) 
ELSE 
EXT(N) = XOEXTR(N) 
EYT(N) = SLOPECN)*(XOEXTR(N) - XOCN)) + YOCN) 
END IF 


* COMPUTATION OF VALUES OF THE VARIABLES NEEDED TO DO LINEAR REGRESSION 
SeenetHe TARGET FIXES TO DETERMINE ESTIMATED COURSE 

LLL = LLL + 1 

SUMX = SUMX + EXT(N) 

SUZ SUR page ICN send 

SUMY = SUbiv = EY RN} 

SUM 2e=  SUMY 2 ey CEYLON) 2 

SUM = SUMKY + EXTCN)itE YIN) 


7.33 CONTINUE 
XMEAN = SUMX/LLL 
YMEAN = SUMY/LLL 


pee tiATED SLOPE (ESLOPE) AND Y INTERCEPT CEYINT) OF TARGET TRACK 
ESLOPE = (SUMXY - SUMX*YMEAN)/(SUMX2 - SUMX*XMEAN) 
EYINT = YMEAN - ESLOPE*XMEAN 


PeoolIMATTON OF COURSE,SPEED AND RANGE OF THE TARGET 
BORER = "(Pi7 tio = ATANCESLOPE)) *DEGPI/P1 


Pee en ee ENG eo OB ShOPEs EC 1) 

PP ee vere 2s tag! ) aoe y een Orr iA 1¢ 27° TE LC I ))) 

N22) 9b Cl 22 Cle Bok PE 32) 

YEP = EY 2°71 ECt) = BRE / Ci spr oO 2 ) 

N22 = 122" EVINT)/ESLOPE 

XFF = (YFF - EYINT)/ESLOPE 

Poke = sORICOXEE = X22) 0250 YOR = Zo)? )/ 
& DRL eine) 

ERANTL = SQRT(CCXFF - XOC2*TELCI)))**2 + 
& CYFER = YOC 247 ECT ee) 

RANGEF = SORTCCSXI( 2*1EECI) = XGC2*TELCI)))**2 + 

& Core z Tir Gl) Saez nl jr) 72 ) 


* INITIAL COMPUTATIONS TO DETERMINE THE PRECISION OF THE PROCEDURE 
CC =CC + 1 


* PRECISION FOR ESTIMATED TARGET RANGE 
PER = (CERANTL - RANGEF)*100. DO) /RANGEF 
Pe CABS CEERI t he 2O a DO iHEN 
VEINTE = VEINTE + 1 
ieCAno (PER 2h, 10,00 THEN 
DIEZ = DIEZ + 1 


IFC ABS PER) Siti en) O)) Bien 
CINCO = CINCO + 1 
END IF 
ENDS 
EGSE 
END IF 


* PRECISION FOR ESTIMATED TARGET SSEEED 
PER1 = (CESPETL* 18000 - Ssack) )*100. 00) 7 ssl 
IFC ABS( PERT) 3 GE 920200) THEN 
VEINT) = VEINt ee 
IFCABSC PER] ) ob vO eo TREN 
DIEZ] = DIEZ 1 
IF CABSCPERI) 2 ee O THEN 
CINCO1 = CINCO1 + 1 
Ea ae 
END er 
ELSE 
ENGR 


* PRECISION FOR ESTIMATED TARGET COURSE 
PER2Z = ECTLR -3se0( J) DEGer en 
IFC ABS(PER2) [Ee 2020p HEN 
VEINT2 = VEINT2 + 1 
TPCABS(CPER2Z) . LE 10.00 tar 
DIEZ2°= DIEZZ ee 
IFC ABSCPER2) [EE o. 2 Op tee 
CINCO2 = CINCO2 + 1 
ptr 1g 
ENE 
ELSE 
ENTREE 


* QUTPUT INDICATING SIMULATED AND ESTIMATED TARGET TRACK PARAMETERS 
WRITEC1, (€13,2X,F9.1,1X,F1201, Zr 5 


& ' )M,RANGEF , ERANTL, SCT(J)**DEGPI/PI,ECTLR,SST(K) , 
& Eoreii loca he 

US CONTINUE 

125 CONTINUE 

850 CONTINUE 

900 CONTINUE 


* FINAL COMPUTATIONS TO DETERMINE THE PRECISION OF THE PROCEDURE. 

R5 = (DBLE( CINCO) /DBLEC CC) )*100 

R10 (DBLE( DIEZ) /DBLEC CC) )*100 

R20 (DBLE( VEINTE ) /DBLE( CC) )*100 
S5 = (DBLECCINCO1)/DBLE( CC) )*100 
510 (DELECDIEGA 7 DELECGE > 100 
S520 (DBLE( VEINT1)/DBLE(CC) )*100 
C5 = (DBLEC(CINCO2)/DBLE(CC) )*100 
C10 (DBLECDIEZ2)/DBLEC CC) )*100 
C20 (DBLE( VEINT2)/DBLE( CC) )*100 


** QUTPUT THAT SHOWS THE PRECISION OF THE PROCEDURE 
WRITE(1,'(C11,1X,F7-1, 15 Foon or 5. 1) 
& M, SO(CL), (SD/PI)“DEGPI, 1YChHMM j= RS, R10.) Reo, 


LA) 
Jes 


& 


Sees On S205 Co ClO 020 


950 CONTINUE 
1100 CONTINUE 


1999 CONTI 

2000 CONTIN 
STOP 
END 


NUE 
UE 


Tedesdedevededcdedeteds ded tededs deve devededevedk ie tesdededaedeveddedetedcdvdedededvdedededvvevedcdvdescakvdkedevedededesvsede deve vedede dex 


** THIS SUBROUTINE GENERATES ONE NORMAL (0,1) NUMBER BY USING THE BOX * 
* MULLER METHOD. THIS NUMBER IS USED BY THE MAIN PROGRAM TO CREATE * 
* THE ERRORS IN THE BEARINGS WHICH ARE ASSUMED NORMAL (0, VARERB) * 
MIeaede ves dene vods veve deve Hae seve He Fe Ved Hee Fe BeBe Fede Fete Tee Te Tete Fes ee Ne TET ETE TCT TTT KIEL TTT RTC TESTER TERRE 

SUBROUTINE NORRN(DSEED, BEERR) 

Een PRECISION UC2), D31M1, DSEED, TWO, Pl, BEERK 

DATA D31M1 /2147483647.D0/, TWO/2.D0/, PI/3. 141592654D0/ 


* LINEAR CONGRUENTIAL METHOD TO OBTAIN A UNIFORM(0,1) 


We 5 I 


=] ,2 


DSEED = DMOD( 16807. DO*DSEED ,D31M1) 


UCI) = DSEED/D31M1 


5 CONTINUE 


* BOX AND MULLER METHOD TO OBTAIN A NORMAL(CO, 1) 
BEERR = (SQRT( - TWO*LOG(UC1))))*COSCTWO*PI*U( 2) ) 


RETURN 
END 


FOSS ella Fak bal ite oe de ee ee ee i i oo i oe a oie el i ie Oe ie aio i ba i bio i ok 
* SUBROUTINE TO COMPUTE THE PARAMETERS OF THE ARCTANGENT MODEL (BO, ¥ 
* TO, AND VD) BY USING FAIRED BEARINGS FROM THE QUADRATIC MODEL. % 
fers. THIS CASE IS USEFUL ONLY FOR TEN MINUTE LEGS. * 


% 


‘ 


DS SISSIES at ah iat rh al letfetel ce |e} ablatyich airtae bar rp bs vee ci sritgr bars tala ira soe) hard ae) oe ke Vor i eke oie ae ae ok rick ak ge he ae | se oe eee eee) eo oe hae) ae a hs 
SUBROUTINE TNGENT( FBE,IN,M,BO,TO,VD) 
PARAMETER(PI = 3. 141592654D0) 
Pence eneolslON BO, £0, VU, AY, YZ, ZA, Tl, 12, T3, FBECIN), 
epee of, FG, GH, SUMTO, SUMVD, SUMBO 


INTEGER 
SUMTO = 0. DO 
SUMVD = 0. DO 
SUMBO = 0. DO 
COUNT = 0 
iran. EQ. 1)THEN 
JK = 49 
KJ = 7/8 
ELSE 
JK = 1 
KJ = 30 
END IF 
POeer iil = JK,(KJ-14) 
AB = III 
e = itl + 7 
Cbe= iil + 14 
T1 = DBLECAB) 
T2 = DBLECBC) 
T3 = DBLE(CD) 
Pie— ©2=<1T3 
Peels 11 


ab epee CD. M, 


IN, COUNT, JK, KJ 


[Nn 
tsAy 


PG =o =) 12 
GH S12 2 Til 
XY = TAN(FBE(BC) - FBE(AB)) 
YZ = TAN(FBE(CD) - FBE(AB)) 
ZA = TAN(FBE(CD) - FBE(BC)) 
TO = (XY*YZ*T1*(FG**2) - ZA*XY*T2*(EF**2) + ZA*YZ*T3*(GH**2))/ 
(DE*XY*(ZA*EF - YZ*FG) - GH*ZA*(XY*EF - YZ*GH)) 
VD = (XY*YZ*DE)/(XY*EF*(T2 - TO) + YZ*GH*(TO - T3)) 
BO = FBE(AB) - ATAN(VD*(T1 - TO)) 
IF(BO .GT. 2.DO*PI)BO = BO - 2.DO*PI 
IF(BO .GT. 3.DO*PI/2.D0)BO = BO - 2.DO*PI 
SUMTO = SUMTO + TO 
SUMVD = SUMVD + VD 
SUMBO = SUMBO + BO 
COUNT = COUNT + 1 
CONTINUE 
TO = SUMTO/DBLE( COUNT) 
VD = SUMVD/DBLE( COUNT) 
BO = SUMBO/DBLE( COUNT) 
RETURN 
END 


("Nn 
ON 


10. 
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